Long-range dynamics of magnetic impurities coupled to a two-dimensional Heisenberg 

ant iferromagnet 
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We consider a two-dimensional Heisenberg antiferromagnet on a square lattice with weakly coupled 
impurities, i.e. additional spins interacting with the host magnet by a small dimensionless coupling 
constant g <gi 1. Using linear spin- wave theory, we find that the magnetization disturbance at 
distance r from a single impurity behaves as 5S^ ~ g/r for 1 <^ r <^ 1/g and as 5S^ ~ 1/(5'*'^) 
for r ^ Surprisingly the magnetization disturbance is inversely proportional to the coupling 

constant! The interaction between two impurities separated by a distance r is 5e oc g"^ /r for 1 ^ 
r <C l/g and 5e <xl/r^ for r 3> l/g. For large distances, the interaction is therefore universal and 
independent of the coupling constant. We have also found that the frequency of Rabi oscillations 
between two impurities is logarithmically enhanced compared to the decay width cjRabi tx g^ ln{l/ gr) 
at 1 <^ r ^ l/g. This leads to a logarithmic enhancement for NMR and EPR line broadening. All 
these astonishing results are due to the gapless spectrum of magnetic excitations in the quantum 
antiferromagnet . 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.30.Hx 



I. INTRODUCTION 

The interplay of magnetic impurities with strongly cor- 
related electron systems has attracted considerable at- 
tention over the past decade. The discovery of high- 
temperature superconductors stimulated studies of impu- 
rities in two-dimensional Mott-insulators with long-range 
antiferromagnetic orderi^, where the copper-oxide par- 
ent compounds are driven to superconductivity by dop- 
ing with holes or electrons. At low doping, the dopants 
are localized and it is therefore insightful to study the 
limit of isolated static holes, which have been realized 
experimentalljiSiii^iSil and extensively studied theoreti- 
calljiSiSiiSiiiiiSiiii^. There is a considerable body of work 
devoted to impurity bonds and added spin e-^'^d^d^d'' , a 
generalization of the static hole case, and an unexpected 
behavior of the impurity magnetic susceptibility at low 
temperature has been revealed quite recentlyiSiiSiSfliSiiS^. 
There is also a separate and very interesting Kondo-like 
problem of an impurity in a magnetic system close to an 
0(3) quantum critical pointiS^iSiSiSl. However, in this 
work, we consider impurities (added spins) in a system 
with long-range antiferromagnetic order. The very un- 
usual behavior we find, closely related to that observed 
for the magnetic susceptibility and other quantities in 
the paper c^d^i^^i^-'^i^^i^^ , is due to gapless Goldstone spin- 
wave excitations. 

The rest of the paper is organized as follows. In section 
^we formulate the model, derive the spin- wave vertices 
and introduce an effective Hamiltonian describing the in- 
teraction between impurities, which is explicitly calcu- 
lated in sections IIIII an IIVI as a function of the distance 
between impurities. A similar problem has been consid- 
ered earlier in Ref. 8, but only the case of small separa- 
tions has been addressed, here we focus on the long-range 
behavior. The magnetization disturbance in the host an- 
tiferromagnet induced by a single impurity is then cal- 



culated in section and finally section VI presents our 
conclusions. 



II. MODEL 

We consider one and two spin — impurities coupled 
to an isotropic Heisenberg antiferromagnet (J > 0) on a 
square lattice with lattice spacing a = 1. The impurities 
are connected to the origin and site r of the antiferro- 
magnet. The Hamiltonian of this system reads 



H 



Sj + j' ((Ti • So + T2 ■ Sr) 



(1) 



Here denote nearest-neighbors, Si is a spin-S* oper- 
ator at site and are spin-i operators describing the 
impurities, which are either both ferromagnetic (J' < 0) 
or antiferromagnetic (J' > 0). After integration over 
quantum fluctuations of the system, one obtains the ef- 
fective Hamiltonian for the interaction between impuri- 
ties 



H, 



eff 



e(r) ala^ 



(A/(r) 



H.c) . (2) 



The calculation of this effective Hamiltonian is the goal of 
the present work. The diagonal interaction term e (r), as 
well as the off-diagonal term M (r) are different, depend- 
ing on whether the impurities are coupled to the same or 
different sublattices, see Fig. QJa) and (b) respectively. 
Without loss of generality, we place the origin on the up 
sublattice (sublattice "a" ) and indicate the sublattice the 
second impurity is coupled to by the corresponding let- 
ter, see Figs.Q] In this way, "a" ("b") refers to impurities 
coupled to the same (different) sublattice. We study the 
weak coupling regime | J'| <C J, so that the dimensionless 
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FIG. 1: Schematic picture of the antiferromagnetic host with 
two antiferromagnetic impurities coupled to the origin and 
site r on the same sublattice (a) and on different sublattices 
(b). 



TABLE I: Spin-wave vertices for an impurity coupled to site 
r on sublattice "a" or "b". 

Symbol Operator Factor for r on "a" or r on "b" 
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coupling constant is small 

g = |J'|/(2%/2J) < 1 



(3) 



□ 



J' 



J\l jfV^e^'^-^ 



J' 
J' 



J V N 



To account for quantum fluctuations, we treat the host 
antiferromagnet in the linear spin-wave approximation 
where excitations are described by operators and 
creating spin- waves with — ~1 and = +1 respec- 
tively, see Ref. Q for a review. In this approximation, 
the Hamiltonian can be decomposed as H — Hq + V , 
where 

CJq (a^Qfq + /3q/3q) + J' S {fff ± (T2 

q 



V 



H 



(a) 



with 



(4) 



Eq — iJSujci- The interaction Hamiltonians © and © 
generate one and two spin-wave vertices summarized in 
Table D 

Let us return to the effective Hamiltonian . The first 
two terms are obvious and do not require calculations. 
They simply generate three energy levels, Eq = — J'(S'^), 
El ^ = 0, and E2 = J'iS""). To be specific, let 
us consider the case shown in Fig. ^a) , with impurities 
coupled to the same sublattice, then 



p,q 

+ ^'/l (^^.+ Ee"' '"'K"q + «q/5q)+h.c.j - (5) 



|0) = 

|1> = 

|1> = 

|2> = 



i,T) , 
T,i> , 
T,T> . 



(7) 



H 



ib) 



p.q 

+ ^'/f('^.+ E^"''"'Mq + «q<)+ll.C.^ . (6) 

Here Eq is the ground state energy of the antiferromag- 
netic host and N the number of sites. The upper (lower) 
sign in the ± expression refers to the situation where the 
second impurity is coupled to a spin on sublattice "a" 
(sublattice "b" ) . This convention is used throughout the 
whole paper. The Bogoliubov parameters Uq and Vq are 
given by 



Only the ground state |0) is the true stationary quan- 
tum state. The states |1), |1), and |2) decay to the 
ground state with emission of spin-waves. Using the 
Fermi golden rule and the decay matrix elements pre- 
sented in Tab.m one finds the following widths of excited 
states with respect to the emission of magnons 



F2 = 2Fi . 



(8) 




with u)q = — and 7q — ^ (cosg^r + cosgy), see 
Ref. Q. In this notation, the spin-wave dispersion is 



The three-level system is well defined, since Ei — Eq 
Fi. But the diagonal and off-diagonal interaction ener- 
gies e (r) and M (r) in the effective Hamiltonian (O have 
limited meaning because of the finite lifetime. Only in the 
ground state |0), the diagonal interaction energy is well 
defined. In the spin flip states and |2) the diago- 

nal interaction energy does not make much sense, because 
we will see that it is always much smaller than the corre- 
sponding decay width. However, there is a regime where 
the off-diagonal interaction M (r) is larger than the de- 
cay width and hence leads to Rabi oscillations between 
states |1) and |T). 



3 





t 


, i 


i 






I 




■^^^ 




* 


t 


t 


T 






i 


— 1 


i 






t 


t 


i 








/ 

^ / 
~- 

/ -~ 

-^^J 










i 1 




i 


i 







i 


\ 








i 


I 


/ 


t t 


X 


i 




i 


— 


i 






i 


\ 

\ 







; 1 \ i 



t Vii^ t 



FIG. 3: Diagrams describing the ^-corrections for antiferro- 
magnetic impurities coupled to sites on different sublattices 
(case "b"). The notation is the same as in Fig.\^ 



FIG. 2; Diagrams describing the ^-corrections for antiferro- 
magnetic impurities coupled to sites on the same sublattice 
(case "a"). The spin-wave vertices are summarized in Tah.\^ 
dashed and dotted lines represent a- and (5 -spin-waves respec- 
tively. 



III. DIAGONAL INTERACTION e(r) BETWEEN 
IMPURITIES 



In the leading order of the ^-expansion, the interaction 
energy e (r) arises in second, third and fourth orders of 
perturbation theory, describing the exchange of two spin- 
waves between impurities. The corresponding contribu- 
tions in usual Rayleigh-Schrodinger perturbation theory 



5^2 

6tA 



E 



E 
E 



{mn){n\V\Q) ^ 

(eo - En) 

{Q\V\n){n\V\m){m\V\Q) 



(eo - en)(eo - Cm) 

{Q\V\n){n\V\m){m\V\k){k\V\Q) 
(eo - e„)(eo - e„i)(eo - e^) 



(9) 



where V is the perturbation (0J. In general, the expres- 
sions for third and fourth order energy corrections are 
more complex than those in The complication is 

due to contributions similar to lO, but in which inter- 
mediate states m,n, k coincide with the initial state |0). 
Fortunately, the perturbation does not allow such in- 
termediate states and therefore Eqs. are valid. The 
matrix elements (n|I/|r7i) are given in Tableland it is 



convenient to represent the corrections © by diagrams 
where each vertex corresponds to some particular ma- 
trix element (ri|V^|m). The diagrams describing the ^- 
corrections are shown in Figs. [21 and |31 for antiferromag- 
netic impurities (J' > 0) coupled to the same and dif- 
ferent sublattices respectively. Ferromagnetic impurities 
( J' < 0) generate similar diagrams with interchanged a- 
and /?-spin- waves. The first diagram in Fig.|21for instance 
represents the expression 

-2J'^^2y2^2g-,(p-q).r 

iV2 (I J'l S -f AJSojp) US (wp + wq) (I J|' S + AJSojp) ' 

with a factor 2 taking into account the similar process 
with exchanged impurities. 

Summing these contributions, we obtain the interac- 
tion energy in the leading -^-order 



5e (r) 



Ar 



o2(p-q)-r 



Br. 



(10) 



with 



("p ('^p ^ ^q 

-|-8Jj'WqWptJq (Wp (Wp + Wq) - 2WptJq) 



-2MpUqWpWqWpWq 



(a) ^ 

p,q 



p,q 



(Wp + Wq) (I J'l + 4Jc^p)' (I J'l + 4JWq)' 

(a;p+Wq)(|J'|-|-4Ju;p)(|J'|-t-4Jc^q) . 
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The sum extends over momenta in the magnetic Brillouin 
zone. In the case of ferromagnetic impurities (J' < 0), 
the Bogohubov parameters in 1)10(1 have to be inter- 
changed. 

For large distances, r ^ 1, the interaction energy 
comes from small momenta, p, g ~ 1/r <C 1. One can 
therefore approximate the Bogoliubov parameters and 
the dispersion by 



q 

V2 



V2q 



and 



V2q ' 



(11) 



and hence simplify the interactions (|10|l 



5e(") (r) 
(5e('') (r) 



Jg^ 2V2 f {g^ - qp) e*(P-i) '-d2p d'^q 



S (27r) 



ip- 



h?) {g+py (g + q) 



S {2nf J ip + q)i9 + p)ig + q) 



(12) 



Here g is the dimensionless coupling constant (Q- Since 
both integrals are ultraviolet convergent, we extend the 
integration domain to infinity. There is no difference be- 
tween ferro- and antiferromagnetic (J' — ±| impuri- 
ties in this approximation. 

We first consider very large distances r ^ 1/g. In this 
case, momenta in (I12|l are limited hy p, q <^ g and the 
corrections to the ground state energy are equal to 



gi(p-q)-i 

(p + q) 



-(Pp d?q , 



where according to our convention the plus sign corre- 
sponds to 5t^°''^ (r) (same sublattice) and the minus sign 
corresponds to (Je*^^^ (r) (different sublattices). A first in- 
tegration over angles gives two Bessel functions Jo, which 
can be integrated using J- Mpr)M,r)p,dp<i, ^ 

yielding 



5e (r) = ± 



1 



J 



16V27r Sr^ 



[r > 1/5) 



(13) 



Interestingly, in this limit, the interaction between im- 
purities is independent of their coupling to the antiferro- 
magnetic host J' . 

In the case of intermediate distances between im- 
purities, 1 <C r <^ 1/g, a similar calculation using 

00 Mpr)Mqr)dpdq ^ 1^ jg^jg 



Jo 



P+Q 



Seiv) 



J 



2V2TTSr 



(1 < r « 1/g) 



(14) 



If the separation between the impurities is small, r « 1, 
the approximation (|ll|l for the Bogoliubov parameters 
and the dispersion is no longer valid and the integrals 
in (|10|) have to be calculated numerically. In order to 
circumvent finite-size effects, we extrapolate the results 
obtained from lattices with up to N"^ = 100 x 100 sites 



by a polynomial Se (r) = a(r) -|- b{r)/N + c{r)/N'^. Even 
for very small distances, the interaction between addi- 
tional spins is surprisingly well fitted by the intermediate 
distance attractive behavior H14|l . The situation is anal- 
ogous to the magnetization disturbance, see section V, 
for which Fig. [S] provides an illustration of qualitatively 
similar fits. The attractive interaction (|14|) found in the 
weak coupling regime is independent of the sublattice. 
In comparison, in the strong coupling limit \J'\ > J an 
attractive interaction is found for nearest neighbors, but 
the interaction between impurities on next-nearest sites 
is repulsiveSiS. 

In the above calculations, we assumed that the im- 
purities are in the ground state, see Eq. Q. To check 
the kinematic structure of the diagonal interaction in the 
effective Hamiltonian ^ one has to perform similar cal- 
culations for states |1), |1), and |2), see ((TJ. In the case 
where both impurities are flipped with respect to their 
ground state configuration, the interaction is the same 
as in the ground state, since the situation corresponds to 
ferromagnetic coupling with J' > 0. For an excited state 
with only one flipped impurity, a calculation analogous 
to (|10|l shows that the interaction has opposite sign and 
hence justifies the kinematic structure of the crjul term 
in |2Il. We emphasize (see also end of section InJ, that 
because of the finite lifetime, the diagonal interaction in 
the spin-flipped states has limited meaning. 



IV. OFF-DIAGONAL INTERACTION M (r) AND 
RABI OSCILLATIONS BETWEEN IMPURITIES 

If a magnetic impurity is flipped with respect to the 
ground state configuration, e.g. in a nuclear magnetic 
resonance (NMR) or electron paramagnetic resonance 
(EPR) experiment, then the excited state has a finite 
lifetime due to emission of magnons, see first of Eqs. JHJ. 
However, there is another mechanism of spin relaxation, 
due to presence of distant similar magnetic impurities: 
a spin-wave exchange leads to Rabi oscillations between 
two impurities. In the effective Hamiltonian ^ this pro- 
cess is described by the off-diagonal term M (r) . Di- 
agrams for M (r) are shown in Fig. 01 Rabi oscillations 
can only be observed between impurities coupled to spins 
on the same sublattice, so 



M^") ( 
M^') (r 



= 



The oscillation frequency of this two-level system is pro- 
portional to the real part of the mixing matrix element 
M = (t, il H li, t), since the probability P (t) to find the 
system in a state with flipped impurities after the time t 
is equal to P (t) = sin^ (Re M t). Using Tab. U to evalu- 
ate the diagrams shown in Fig. ^ we find 



J" 



4 JWq -I- iS 



J' 4 Jlj„ 



5 



J 



t 1_ i 

I y\ t 

FIG. 4: Diagrams describing Rabi oscillations: The system 
is brought into a state where one spin is flipped, for exam- 
ple during an NMR experiment and oscillates between excited 
states by exchanging spin-waves. 



For large distances r ^ 1 we use the approximate Bo- 
goliubov parameters The mixing element becomes 



M^'^) (r 



Ja (qr) 



2q 



q- - 



+ iT^S {q- g)] dq 



Jg^ {Yo {gr) - iJo (gr)) 



where lo is the Neumann function. For very large dis- 
tances between impurities, r 3> l/g, the real part of the 
mixing matrix element is comparable to the imaginary 
one and both are much smaller than the width (jSJ. In 
this case, there are no Rabi oscillations. However, in the 
intermediate regime, 1 ^ r ^ the real part 



Re M'-"^ (r) = Jg^-ln 



gr 



(15) 



is logarithmically enhanced compared to the imaginary 
part and compared to the spin- wave width ©. Thus, in 
this regime, Rabi oscillations between impurities are well 
pronounced and this mechanism gives the main contri- 
bution to the effective width of the magnetic resonance 
line 



4^2 



(16) 



V. 



MAGNETIZATION CLOUD AROUND AN 
IMPURITY 



An interesting question is how an additional spin in- 
fluences the magnetic order in the host antiferromagnet, 
described by the staggered magnetization. Let us con- 
sider an impurity ai at the origin and a local magnetic 
field h on site r. The Hamiltonian reads 



H 



S, -f J'ai ■ So + hS^ 



(17) 



where 5E is the part of the energy dependent on J'. Now 
we can consider the same framework as for the calcula- 
tion of the interaction between impurities. Using Eq. 
(|18|l together with the Rayleigh-Schrodinger perturba- 
tion theory, one can see that 

SS^ = T2Se (r) , 

where 6e is the energy correction described by those di- 
agrams in Figs. 12 and 121 without spin flip of the second 
impurity. These diagrams give the following explicit ex- 
pressions for the magnetization variation 



=«(p-q)-r 



p>q 



D 



(19) 



with 
^vn - 

^P.q = 



2j'Mq'i;pWq - J'UpUq (Wp + CJp) + SJUq-WpWpCJq 
-SJUpUqVpWqCJpLJq 

(C^p + C^q) (I J'l + 4Jc^p) (I J'l + 4Jc^q) 



assuming J' > 0. The expressions for the ferromagnetic 
case (J' < 0) are obtained by a change of the global sign 
and interchanged Bogoliubov parameters in (|19|) . Using 
the approximations (|ll|l for r 3> 1 , we find that the mag- 
netization disturbance is equal to 



8St 



±- 



1 J 1 



= ±- 



1 1 



'8%/27r |J'| S'r3 ^■i2'KgSr'^ 
in the very large distance limit, and 



(r > 1/5) 



(20) 



8Sl 



± 



1 U'l 1 



± 



5 



SVStt J Sr Att Sr 



(1 < r< 1/5) 



(21) 

for the intermediate region. In the case of small separa- 
tion from the impurity, one has to calculate the integrals 
in 1)19(1 numerically. Using the same finite-size extrapo- 
lation scheme as in section HTll wc find that the variation 
of the magnetization in the vicinity of the impurity is 
described by Eq. (|^ down to r = 1. Fig. displays the 
variation of the magnetization calculated numerically at 
J' = —0.01. The results are very close to the analytical 
expression (|21f) . In agreement with Ref . H an added spin 
always enhances the Neel order in the host magnet, inde- 
pendent of the sign of the exchange coupling. In contrast 
to a vacancy, this enhancement is not limited to nearest 
neighbor sites^iiij but extends over the whole magnet. 
It is also interesting to compare our result for an added 
spin to an in-plane impurity considered in Ref. Hi. If the 
impurity is placed inside the host, it weakens the Neel 
order of the surrounding spins, but the magnetization 
disturbance also decreases as 1/r^ for r ^ l^**. 



and the variation of the expectation value of the local 
magnetization dS^ in the ground state is given by 



6S^ 



dSE 
~dh 



VI. CONCLUSION 



To conclude, we have studied the long-range dynam- 

(18) 

ics of one and two spin-i impurities in a two-dimensional 
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Sublattice "a" 

Fit: 5.63- lO'V (r>2) 

Sublattice "b" 

Fit: -5.94- 10"V (r>l) 



FIG. 5: The points show the results of a numerical calculation 
of the variation of the host magnetization (S = \) at site r 
induced by a ferromagnetic impurity coupled to the origin by 
J' — —0.01. The curves show fits of these points according to 
SSr ~ F/r, where F is the fitting coeffictent. Fitted values of 
F (see inset) are very close to the value F — 5.62 • 10~^ that 
follows from Eq. 



Heisenberg antiferromagnet with on site spin-S" treated in 
the Hnear spin-wave approximation. The impurities are 



assumed to be weakly coupled to the host magnet by a 
small dimensionless coupling constant g. A systematical 
treatment of the corrections contributing to the leading 
order of the ^-expansion leads to non-trivial long-range 
dynamics. The interaction between two impurities can 
be separated into two regimes: For very large separa- 
tions (r ^ 1/5) it is universal (independent of J') and 
decreases as \/r^. The interaction is repulsive (attrac- 
tive) for impurities coupled to the same (different) sub- 
lattices. In an intermediate region 1 <C r <C the 
interaction decreases only as 1 /r and is attractive, inde- 
pendent of the sign of the exchange couplings. It is shown 
that Rabi oscillations between impurities coupled to spins 
on the same sublattice are possible and well pronounced 
in the intermediate regime. The effective Hamiltonian 
decsribing the interaction in terms of the impurity spins 
is derived. It exhibits an xyz anisotropy which leads to 
NMR and EPR line broadening. The magnetization dis- 
turbance in the host magnet induced by a single impurity 
is analyzed in the same framework. It is shown that the 
disturbance exhibits behaviors similar to the interaction 
energy, always enhancing the magnetic order in the anti- 
ferromagnetic host. Numerical results indicate that the 
intermediate regimes can be extended down to r « 1. 
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